ENTROPY PRODUCTION IN GAUSSIAN THERMOSTATS 
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Abstract. We show that an arbitrary Anosov Gaussian thermostat on a surface is 
dissipative unless the external field has a global potential. This result is obtained by 
studying the cohomological equation of more general thermostats using the methods 
in@. 



1. Introduction 

Gaussian thermostats provide interesting models in nonequilibrium statistical me- 
chanics jHl HE E] • Given a closed Riemannian manifold (M, g) and a vector field E 
(the external field) on M, the Gaussian thermostat (or isokinetic dynamics, cf. [T3] ) 
is given by the differential equation 

This equation defines a flow on the unit sphere bundle SM of M which reduces to 
the geodesic flow when E = 0. 

In general, Gaussian thermostats are not volume preserving and the purpose of 
the present paper is to characterize precisely those Anosov Gaussian thermostats in 
2 degrees of freedom which do not preserve any smooth measure. 

When <f) is Anosov and dimM = 2 a result of E. Ghys ^1] ensures that <fi is 
topologically conjugate to the geodesic flow of a metric of constant negative curvature 
and thus cf) is transitive and topologically mixing. For such a flow it is well known 
(cf. [HI Chapter 20]) that there exists a unique Gibbs state p associated with the 
Holder continuous potential — | f log J", where J" is the unstable Jacobian of cf). 
The measure p is characterized by being the maximum of 



h, 



log J? dv 



't 

=o 



where v runs over all 0- invariant Borel probability measures and h u (<f)) is the measure 
theoretic entropy of 4> with respect to v. The unique measure p is called the SRB 
measure of <p. If r is a probability measure which is absolutely continuous with respect 
to the Liouville measure of SM, then p is also the weak limit of ^ (jfr dt as T — >• oo. 
The entropy production of the state p is given by (cf. 



e^(p) := — J div F dp = — Lyapunov exponents 

where F is the infinitesimal generator of <fi and divF is the divergence of F with 
respect to any volume form in SM. 

l 
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Fix a volume form on SM. Any other volume form can be written as /0 for 
some smooth positive function /. If we let LpQ be the Lie derivative of along F, 
then 

L F (fQ) = d(i F f&) = F(f)Q + fL F Q = F(f)G + fdwFQ. 
Hence if divF denotes the divergence of F with respect to /0 we have 

(2) divF = F(log /) + divF. 

In other words the two divergences are flow cohomologous (and thus is well defined 
for any 0-invariant measure). 

Ruelle [20! has shown that e^(p) > with equality if and only if p is also the 
SRB measure of the flow </>_ t . If p is an SRB measure for both <p t and <fr-t then 
the theory of Gibbs states for Anosov flows (cf. fj| Proposition 20.3.10]) implies 
that — ^ 1 1=0 log J i and ^| log J t s are cohomologous (and the coboundary is the 
derivative along the flow of a Holder continuous function). It follows that <fi preserves 
an absolutely continuous invariant measure with positive continuous density (and 
this measure would have to be p). An application of the smooth Livsic theorem ^SJ 
Corollary 2.1] shows that <ft preserves an absolutely continuous invariant measure with 
positive continuous density if and only if preserves a smooth volume form. Using 
(0) we see that e^(p) = if and only if divF is a flow coboundary and we can take 
divF with respect to any volume form. 

Let 9 be the 1-form dual to E, i.e., 9 x [y) = (E(x),v). An easy calculation (see 
Lemma I3.2J1 shows that if we consider in SM the volume form determined by the 
canonical contact 1-form, then divF(x,v) = —0 x (v). Thus e^(p) = if and only if 
there is a smooth solution u to the cohomological equation 

(3) F{u) = e. 

We will show as a consequence of a more general result to be stated below that if 
dimM = 2 then Q holds if and only if 6 is an exact form, i.e. if and only if E has a 
global potential. Thus we obtain: 

Theorem A. An Anosov Gaussian thermostat on a closed surface has zero entropy 
production if and only if the external field E has a global potential. 

A system with e^(p) > is referred to as dissipative. Dissipative Gaussian ther- 
mostats provide a large class of examples to which one can apply the Fluctuation 
Theorem of G. Gallavotti and E.G.D. Cohen |HJ 03 (extended to Anosov flows by 
G. Gentile |10j ) and this theorem is perhaps one of the main motivations for determin- 
ing precisely which thermostats are dissipative. Observe that Gaussian thermostats 
are reversible in the sense that the flip (x, v) \— > (x, —v) conjugates 4>t with (just 
as in the case of geodesic flows). We recall that the chaotic hypothesis of Gallavotti 
and Cohen asserts that for systems out of equilibrium, physically correct macroscopic 
results will be obtained by assuming that the microscopic dynamics is uniformly hy- 
perbolic. 

In |22J, M. Wojtkowski proved Theorem A assuming that E has a local potential 
(i.e. 9 is closed) and in [T|, F. Bonetto, G. Gentile and V. Mastropietro proved the 
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theorem for the case of a metric of constant negative curvature and 9 a harmonic 
1-form. We emphasize that we do not make any assumptions on g or E except that 
the underlying isokinetic dynamics is Anosov. Conditions under which the Anosov 
property holds have been given in EHj • 

We now explain for which Anosov systems we can understand the cohomological 
equation (j3J) completely. 

Let M be a closed manifold endowed with a Riemannian metric g. We consider 
a generalized isokinetic thermostat. This consists of a semibasic vector field E(x,v), 
that is, a smooth map TM 3 (x,v) i— > E(x,v) G TM such that E(x,v) G T X M for 
all (x, v) G TM. As before the equation 

dt l7 ' 7J hf T 

defines a flow on the unit sphere bundle SM. These generalized thermostats are 
no longer reversible unless E(x,v) = E(x, —v). 

Suppose now that M is a closed oriented surface. We can write 

E(x, v) = k(x, v)v + X(x, v)iv 

where i indicates rotation by 7r/2 according to the orientation of the surface and k 
and A are smooth functions. The evolution of the thermostat on SM can now be 
written as 

(4) -_ = A(7, 7) ry. 

If A does not depend on v, then is the magnetic flow associated with the magnetic 
field Af2 a , where Q a is the area form of M. Of course, magnetic flows are Hamiltonian. 
If A depends linearly on v, we obtain the Gaussian thermostat ((TJ. 
Let 7r : SM — > M be the canonical projection. 

Theorem B. Let M be a closed oriented surface and consider a generalized isokinetic 
thermostat Q). Suppose the flow <fi is Anosov and let F be the vector field generating 
cf). Let h G C°°(M) and let 6 be a smooth 1-form on M. Then the cohomological 
equation 

F(u) = h07T + 6 

has a solution u G C° C '{SM) if and only if h = and 9 is exact. 

Note that by the smooth Livsic theorem ^2] saying that h o n + 9 = F(u) is 
equivalent to saying that h o n + 9 has zero integral over every closed orbit of (p. 

Theorem B was proved in [H] for the case of magnetic flows (i.e. A depends only 
on x). It was surprising for us that the theorem also holds for systems that do not 
preserve a smooth measure. The proof is also based on establishing a Pestov identity 
as in |2J m f° r geodesic flows, but some unexpected cancellations take place producing 
in the end formulas which are just what one needs to prove the theorem. Earlier 
proofs of Theorem B for some geodesic and magnetic flows using Fourier analysis can 
be found in P2||IB|. 
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Finally we note that Theorem A also holds if we allow magnetic forces. Indeed 
Theorem B holds for a generalized thermostat and divF = — 9 even when we have a 
magnetic field present. The extension of Theorem A to isoenergetic thermostats (i.e. 
in the presence of potential forces) is discussed in Remark |5. 11 

Acknowledgements: The first author thanks the Department of Pure Mathematics 
and Mathematical Statistics at the University of Cambridge and Trinity College for 
hospitality and financial support while this work was in progress. 



2. Preliminaries 

Let M be a closed oriented surface, SM the unit sphere bundle and 7r : SM — > M 
the canonical projection. The latter is in fact a principal S' 1 -fibration and we let V 
be the infinitesimal generator of the action of S 1 . 

Given a unit vector v G T X M, we will denote by iv the unique unit vector orthogonal 
to v such that {v, iv} is an oriented basis of T X M. There are two basic 1-forms a and 
P on SM which are defined by the formulas: 

The form a is the canonical contact form of SM whose Reeb vector field is the 
geodesic vector field X. The volume form a A da gives rise to the Liouville measure 
d\i of SM. 

A basic theorem in 2-dimensional Riemannian geometry asserts that there exists 
a unique 1-form if) on SM (the connection form) such that i/j(V) = 1 and 

(5) da — ip A (3 

(6) dp = —ip A a 

(7) dil) = -(Kon)aAP 

where K is the Gaussian curvature of M. In fact, the form ip is given by 




where Z : (— e, e) — > SM is any curve with Z(0) = (x, v) and Z(0) = £ and is the 
covariant derivative of Z along the curve n o Z. 

For later use it is convenient to introduce the vector field H uniquely defined by the 
conditions P(H) = 1 and a(H) = ip(H) = 0. The vector fields X, H and V are dual 
to a, P and ip and as a consequence of (J3H7J) they satisfy the commutation relations 

(8) [V,X]=H, [V,H) = -X, [X,H] = KV. 



Equations (JSHZJ also imply that the vector fields X, H and V preserve the volume 
form a A da and hence the Liouville measure. 
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3. AN INTEGRAL IDENTITY 

Henceforth (M, g) is a closed oriented surface and X, H, and V are the same vector 
fields on SM as in the previous section. 

Let A be the smooth function on SM given by (J3J), and let 

F = X + XV 

be the generating vector field of the generalized thermostat. 
From (jHJ) we obtain: 

[V, F] = H + V(X)V, [V, H] = -F + XV, [F, H] = -XF + (K - H(X) + X 2 )V. 

Lemma 3.1 (The Pestov identity). For every smooth function u : SM ->Rwe have 

2Hu ■ VFu = {Fuf + (Huf - (K - H(X) + X 2 ){Vuf 

+ F(Hu ■ Vu) + V{X)Hu ■ Vu - H(Fu ■ Vu) + V(Fu ■ Hu). 

Proof. Using the commutation formulas, we deduce: 

2Hu ■ VFu -V(Hu- Fu) 

= Hu ■ VFu - VHu ■ Fu 

= Hu ■ (FVu + [V, F}u) - Fu ■ (HVu + [V, H}u) 

= Hu ■ (FVu + Hu + V(X)Vu) - Fu ■ (HVu - Fu + XVu) 

= (Fu) 2 + (Hu) 2 + (FVu)(Hu) - (HVu)(Fu) - XFu -Vu + Hu- V(X)Vu 

= (Fu) 2 + (Hu) 2 + F(Vu ■ Hu) - H(Vu ■ Fu) - [F, H}u ■ Vu 

- XFu -Vu + Hu- V(X)Vu 

= (Fu) 2 + (Hu) 2 + F(Vu ■ Hu) + V(X)Hu ■ Vu - H(Vu ■ Fu) 

- (K ~H(X) + X 2 )(Vu) 2 

which is equivalent to the Pestov identity. □ 

Now let O := a A da. This volume form generates the Liouville measure d[i. 
Lemma 3.2. We have: 

(9) LpQ = V(X)&; 

(10) L H Q = 0; 

(11) L v Q = 0. 

Proof Note that for any vector field Y, L Y Q = c?(iyO). Since iyQ = —a A j3 — 
—7T*Q a , where Q a is the area form of M, we see that LyQ = 0. Similarly, Lx& = 
L H Q = 0. Finally L F = L X Q + L XV Q = d(i xv Q) = V(X)Q. 

□ 
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Below we will use the following consequence of Stokes theorem. Let N be a closed 
oriented manifold and a volume form. Let X be a vector field on N and / : iV — > R 
a smooth function. Then 

(12) f X(f)& = - [ fL x B. 

J N J N 

Integrating the Pestov identity over SM against the Liouville measure du, and 
using (fTUj) and (fTTj) we obtain: 

2 / Hu-VFudu = / (Fu) 2 du + / (Hu) 2 dfi 

J SM J SM J SM 

(F(Hu ■ Vu) + V(X)Hu ■ Vu) dfi 

SM 

{K-H{X) + X 2 ){Vu) 2 dfi. 

SM 

Using (P2J) and © we get: 

/ (F(Hu ■ Vu) + F(A)iy« • Vu) dfi = 

J SM 

and thus 

(13) 2 / Hu-VFudfi= / {Fu) 2 dfi+ / {Huf d^ 

JSM J SM JSM 

- [ (K -H(X) + X 2 )(Vu) 2 dfi. 

J SM 

We will derive one more integral identity. By the commutation relations, we have 

FVu = VFu -Hu- V(X)Vu. 

Therefore, 

(FVu) 2 = (VFu) 2 + (Hu) 2 + (V(X)) 2 (Vu) 2 
-2VFu ■ Hu - 2VFu ■ V(X)Vu + 2V(X)Vu ■ Hu. 
Thus using again the commutation relations: 

(FVu) 2 = (VFu) 2 + (Hu) 2 + (V(X)) 2 (Vu) 2 

-2VFu ■ Hu - 2FVu ■ V(X)Vu - 2(V(X)) 2 (Vu) 2 . 



Since 

we obtain: 



F(V(X)(Vu) 2 ) = 2V(X)Vu ■ FVu + (Vu) 2 F(V(X)) 

(FVu) 2 = (VFu) 2 + (Hu) 2 - (V(X)) 2 (Vu) 2 
-2VFu ■ Hu - F(V(X)(Vu) 2 ) - (Vu) 2 F(V(X)). 
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Integrating this equation we obtain: 

(14) 2 / Hu-VFudp = [ (VFu) 2 dp + [ (Hu) 2 dfi 

J SM J SM J SM 

- [ {FVufdp- [ F{V{\)){Vu) 2 dp 

J SM JSM 

since by ()12j) and (jSJ) we get: 

/ {F(V(\)(Vuf) + (V(\)) 2 (Vu) 2 } dp: = 0. 

J SM 

Combining (J13)) and (|14j) we arrive at the final integral identity of this section: 
Theorem 3.3. 

(15) I {FVufdp- [ K(Vu) 2 dp = [ {VFufdp- [ {Fuf dp, 

JSM JSM JSM JSM 

where K := K - H(X) + A 2 + F(V(X)). 

Of course this identity holds without any assumption on the underlying dynamics. 
In the next section we will show how to use the Anosov hypothesis to rewrite the left 
hand side of (fT5j) in terms of the stable or unstable bundles. At this point the proof 
differs from the one presented in [3]. We can no longer estimate the left hand side 
of ([l5j) using closed orbits and the non-negative Livsic theorem jTHl EI] since in our 
context the Liouville measure is not necessarily invariant. 

4. Using the Anosov property 

Recall that the Anosov property means that T(SM) splits as T(SM) =M,F®E U ® 
E s in such a way that there are constants C > and < p < 1 < 77 such that for all 
t > we have 

||#-iU«|| < Crf 1 and < C p*. 

The subbundles are then invariant and Holder continuous and have smooth integral 
manifolds, the stable and unstable manifolds, which define a continuous foliation with 
smooth leaves. 

Let us introduce the weak stable and unstable bundles: 

E + = RF © E s , 
E~ = RF © E u . 
Lemma 4.1. For any (x,v) G SM, V(x,v) ^ E ± (x,v). 

Proof. Let A(SM) be the bundle over SM such that at each point (x,v) G SM 
consists of all 2-dimensional subspaces W of T( XtV )SM with F(x, v) G W. 

The map (x, v) 1— > V := RF(x, v) © MV(x, v) is a section of A(SM) and its im- 
age is a codimension one submanifold that we denote by Ay. Similarly the map 
(x,v) 1 — > RF(x,v) © RH(x,v) is a section of A(SM) and its image is a codimension 
one submanifold that we denote by Ah- 
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The flow cj) naturally lifts to a flow <p* acting on A(SM) via its differential. Let F* 
be the infinitesimal generator of <fi* . 
Claim. F* is transversal to Ay. 

To prove the claim we define a function m : A(SM) \ Ah — > K as follows. If 
W G A(SM) \ Ah, then if ^ VF. Thus there exists a unique m = m(W / ) such 
that rriH + V G W . Clearly m is smooth and Ay = m _1 (0) C A(SM) \A H . Fix 
(x, u) G 5M and set m(t) := m((f)^(V(x,v))). By the definition of m, there exist 
functions x(t) and such that 

m(t)H(t) + V(t) = x(t)F{t) + y{t)d(f) t {V). 

Equivalent ly 

m(t)d<j>- t {H{t)) + d4-t(V(t)) = x(t)F + y(t)V. 
Differentiating with respect to t and setting t = (recall that m(0) = 0) we obtain: 

th(0)H + [F, V] = x(0)F + y(0)V. 

But [V, F] = H + V(\)V. Thus m(0) = 1 which proves the Claim. 

From the Claim it follows that Ay determines an oriented codimension one cycle in 
A(SM) and by duality it defines a cohomology class m G J ff 1 (A(S'M), Z). Set E = E ± . 
Given a continuous closed curve a : S 1 — > S'M, the zn<iex of a is i/(a) := (m, [E o a]) 
(i.e. = E*m G i/ 1 (S'M, Z)). The index of a only depends on the homology class of 
a. Since E is 0-invariant, the Claim also ensures that if 7 is any closed orbit of 0, 
then u(j) > 0. 

Recall that according to Ghys we know that <fi is topologically conjugate to 
the geodesic flow of a metric of constant negative curvature. In particular, every 
homology class in Hi(SM, Z) contains a closed orbit of <p. Thus v must vanish. 

If there exists (x, v) G SM for which V(x, v) G E(x, v), then using that every point 
of <p is non- wandering, we can produce exactly as in [T71 Lemma 2.49] a closed curve 
a : S 1 — >• S'M with u(a) > 0. This contradiction shows the lemma. 

□ 

Remark 4.2. The reader will recognize that the index that appears in the proof 
of the lemma reduces to the Maslov index when <p is Hamiltonian. The proof of 
the lemma also follows the presentation in [T71 Chapter 2] of analogous results for 
geodesic flows. 

The lemma implies that there exist unique continuous functions r ± on SM such 
that 

H + r + V G E + , 
H + r-V G E-. 

Note that the Anosov property implies that r + 7^ r~ everywhere. Below we will need 
to use that the functions satisfy a Riccati type equation along the flow. Note that 
are smooth along because E^ are 0- invariant. 
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Lemma 4.3. Let r = r ± . Then 

F(r - V(X)) + r(r - V(X)) + K = 0. 

Proof. Let F = E ± . Fix (x, t>) G SM, flow along and set 

£(t) :=#_ t (#(t) + r(f)V(t)). 

By the definition of r, £(t) G E(x,v) for all i. Differentiating with respect to t and 
setting t = we obtain: 

f(0) = [F,ff] + F(r)y + r[F,V]. 

Using that 

[V, F] = F + y (A)y, [F, H] = -XF + (K- H(X) + X 2 )V 

we have 

f (0) = -\F -rH+ {F(r) + K- H(X) + A 2 - V(X)r} V. 
Replacing H by £ (0) — rU yields: 

6(0) + r£(0) + AF = {r 2 + F(r) + K - H(X) + A 2 - U(A)r} V. 
Since £(0) + r£(0) + AF G F we must have 

r 2 + F(r) + F - F(A) + A 2 - V(A)r = 
which is the desired equation since K = K — H{X) + A 2 + F(V(A)). 

□ 

Here is the main result of this section: 

Theorem 4.4. Let ip : SM — > K 6e a smooth function and suppose <fi is Anosov. 
Then for r = r 

[ {FijjfdjJL- I K^ 2 dfi= [ [F(tP)-rtP + tpV(X)} 2 dfx>0. 

J SM J SM JSM 

Moreover, 

[ [F(ij) -rtP + tpV(X)] 2 dfx = 

J SM 

if and only if ip = 0. 

Proof. Let us expand [F(tp) — rip + -^V^A)] 2 : 

[F(V>) - r</> + </>U(A)] 2 = [F(^)] 2 + ^ 2 r 2 + ^ 2 [U(A)] 2 

- 2F(^)^r + 2F(V')VV r (A) - 2^ 2 rU(A). 

Using that (see Lemma F4.3f) 

F(r - V(A)) + r(r - U(A)) + K = 
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we obtain: 

[F(ip) - rip + ipV(X)} 2 = [F(iP)} 2 - Kip 2 

- F((r - V(X))ip 2 ) + ip 2 [V(X)} 2 - i?rV{\). 

If we integrate the last equality with respect to the Liouville measure \x we obtain as 
desired: 



(Fi/j) 2 dfi- Kip 2 d\i= \ [F(ip) - rip + ipV(X)} 2 dfx 

JSM JSM 

since by (I12|) and Q we have the following cancellation: 

{-F((r - V(X))ip 2 ) + ip 2 [V(X)} 2 - ip 2 rV(X)} dfi = 0. 

■M 

Suppose now 

I [F(ip) -rip + i/jV(X)] 2 dfi = 

JSM 

which implies 

F(ip) -rip + ipV(X) = 
everywhere. Since this holds for r = r^ 1 we deduce: 

(r + — r~)ip = 0. 
But for an Anosov flow r + — r~ ^ everywhere, thus ip = 0. 

5. Proof of Theorem B 



□ 



Let us now prove Theorem B. If Fu = h o n + 9, then it is easy to see that the 
right-hand side of (|15J1 is nonpositive. Indeed, since \x is invariant under v i— > — 1> and 
f — > to we have 

/ ^(u)d/i = and / (^(t>)) 2 rf / u= / (9 x (iv)) 2 d/i. 

JSM JSM JSM 

But V.Fw = x (iv) and thus 

/ (VFu) 2 dfi- [ (Fu) 2 dfi = -[ (hon) 2 d^<0. 

JSM JSM JSM 

Setting ip = Vu, we get 

(16) / {(Fip) 2 — Kip 2 ] dfi < 0. 

JSM 

By Theorem 14.41 this happens if and only if ip = 0. This would give Vu = 0, which 
says that u = f o n where / is a smooth function on M. But in this case, since 
d7T( x ,v){F) = v we have Fu = df x (v). This clearly implies the claim of the theorem. 
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Remark 5.1. Suppose that we include potential forces in our dynamics, that is, we 
consider the isoenergetic thermostat 

(17) ^ = -VW' + E(7 )-^p>7 

on the energy level \\v\ 2 + W(x) = k (we assume that \v\ does not vanish on the 
energy level). Wojtkowski has pointed out [2U Theorem 2.4] that the dynamics of 
(J 17)) reparametrized by arc-length defines a flow on SM which coincides with the 
isokinetic thermostat with external field 

-\/W + E 1 E 
E := g-g = ^V(log(A; - W)) + 



2(k-W) 2 v &v " 2(k-W) 

Since the vanishing of entropy production and the Anosov property are unaltered 
by smooth time changes we conclude applying Theorem A to E that an Anosov 
isoenergetic thermostat has zero entropy production if and only if E/2(k — W) has a 
global potential. 

The question of whether Theorem B extends to higher dimension is more delicate. 
We hope to discuss this topic elsewhere. 
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